Theory of magnetic excitations in iron-based layered superconductors 
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Based on the effective four-band model we analyze the spin response in the normal and super- 
conducting states of the Fe-pnictide superconductors. While the normal state spin excitations are 
dominated by the continuum of the interorbital antiferromagnetic fluctuations and the intraband 
spin density wave fluctuations, the unconventional superconductivity yields different feedback. The 
resonance peak in form of the well-defined spin exciton occurs only for the interband scattering 
at the antiferromagnetic momentum Clj^pj^j for the s± (extended s-wave) superconducting order 
parameter and it disappears rapidly for q < CIafm- The resonance feature is extremely weak for 
the d3,2_y2-wave order parameter due to specific Fermi surface topology of these compounds. The 
essential difference between s±-wave and d^2_y2-wave symmetries for the magnetic excitations can 
be used for experimental determination of the superconducting wave function symmetry. 
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The relation between unconventional superconductiv- 
ity and magnetism is one of the most interesting top- 
ics in the condensed matter physics. In contrast to the 
usual electron-phonon mediated superconductors where 
the paramagnetic spin excitations are suppressed below 
superconducting transition temperature due to the for- 
mation of the Cooper pairs with total spin = 0, in un- 
conventional superconductors, such as layered cuprates 
or heavy fermion superconductors, a bound state (spin 
resonance) with a high intensity forms below T^X^i^. 
The fact that the superconducting gap is changing sign 
at a different parts of the Fermi surface together with 
a presence of the strong electronic correlations yields 
such an enhancement of the spin response^. Most in- 
terestingly, an observation of the resonance peak indi- 
cates not only that Cooper-pairing is unconventional but 
also that the magnetic fluctuations arc most relevant for 
superconductivity^. 

Since the discovery of superconductivity in the quater- 
nary oxypnictides LaFePO^ and LaNiPO^, a new class 
of high-Tc materials with Fe-based layered structure is 
emerging^'^ i^°'^^i^^d3ii4 ^ Although the microscopic na- 
ture of superconductivity in these compounds remains 
unclear at present, certain aspects have been already 
discussedi^'i^^'i^^'20^'22^^^^^. In particular, 
ah initio band structure calculation a^^d^d'^d^dQi^o have 
shown that the conductivity and superconductivity in 
these systems are associated with the Fe-pnictide layer, 
and the electronic density of states (DOS) near the Fermi 
level shows maximum contribution from the Fe-3d or- 
bitals. The resulting Fermi surface consists of two hole 
{h) and two electron (e) pockets. The normal state mag- 
netic spin susceptibility determined from these bandsi22, 
exhibits both small q fluctuations and antiferromag- 
netic commensurate spin density wave (SDW) peaks. 

In this Rapid Communication, using the four-band 
tight-binding model we study theoretically the spin re- 



sponse in the normal and superconducting states of Fe- 
based pnictidc superconductors. We show that the re- 
sulting magnetic fluctuation spectrum calculated within 
random-phase approximation (RPA) consists of two con- 
tributions. The first one is from the antiferromagnetic 
(AFM) spin fluctuations peaked at Qafm = (ti'iTt) aris- 
ing due to the interband scattering. The second contribu- 
tion comes from the intraband scattering and results in 
a broad continuum of the SDW fluctuations with a small 
momenta. We show that the unconventional supercon- 
ductivity yields different feedback on the magnetic exci- 
tation spectrum. The resonance peak in form of the spin 
exciton occurs only for the interband scattering at the 
AFM momentum for the s-t-wave superconducting or- 
der parameter. This peak is confined to the AFM wave 
vector and disappears rapidly away from it. We suggest 
that the superconductivity is most likely s±-wave and is 
driven by the repulsive interaction. 

The Fe ions form a square lattice in the FeAs layer 
of LaFeAsO system, which is interlaced with the second 
square lattice of As ions. Due to the fact that As ions 
sit in the center of each square plaquette of the Fe lattice 
and are displaced above and below the Fe plane, the crys- 
tallographic unit cell contains two Fe and two As ions. 
The band structure calculation o^^'^^i^'^'dSd^i^o show that 
three Fe-3c? states [dxz, dyz, and d^y) give the main con- 
tribution to the density of states close to the Fermi level 
and that these states disperse weakly in the z-direction. 
The resulting Fermi surface consists of two hole (h) pock- 
ets centered around the F = (0,0) point and two electron 
(e) pockets centered around the M = (tt, tt) point of the 
folded Brillouin zone (BZ)^^. Note, the folded BZ corre- 
sponds to the case of two Fe atoms per unit cell, and the 
wave vector (tt, tt) in the folded BZ corresponds to the 
(tt, O) wave vector in the unfolded BZ (related to the case 
of one Fe per unit cell). To model the resulting band 
structure we assume the following single-electron Hamil- 



2 



tonian 
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Qfi, a2, /3i, /32 refer to the band indices, 



are the on-site single-electron energies, t. 



{cos kx + cos ky) + ^2^' ^ cos /cj: cos fcy is the elec- 
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tronic dispersion that yields hole pockets centered around 
the r point, and t^^'^^ = t^^'^^ {cos + cos ky) + 

fcx - - - 



ij^'^^ cos ^ cos ^ is the dispersion that results in 
the electron pockets around the M point. Using 
the abbreviation {e'^,t\,t2) we choose the parameters 
(-0.60,0.30,0.24) and (-0.40,0.20,0.24) for the ai 
and a2 bands, respectively, and (1.70,1.14,0.74) and 
(1.70,1.14,-0.64) for the /3i and (32 bands, correspond- 
ingly (all values are in eV). 

In Figs.[IJa) andj^a) we show the resulting band struc- 
ture and the corresponding Fermi surface for the undoped 
case, X = 0. The Fermi surface consists of the two hole 
(ofi and and two electron (/3i and f32) pockets. The 
(3 bands show much broader bandwidth and are degen- 
erate along X — M direction which is a consequence of 
the hybridization of the underlaying dxz and dyz orbitals 
within the folded BZ. The a bands centered around the 
r point are narrower which also results in the signifi- 
cant contribution to the density of states. The chosen 
band structure reproduces correctly the local-density ap- 
proximation (LDA) Fermi surface topology and the corre- 
sponding values of the Fermi velocities. In particular, we 
have selected the on-site energies and the hopping matrix 
elements assuming the compensated metal at zero dop- 
ing and the filling factor n — A (we further assume that 
there exists another band below the Fermi level which 
is fully occupied and not considered here). Additionally, 
we take into account the details of the electronic disper- 
sions of the bands which form the corresponding hole and 
electron Fermi surface pockets. The comparison between 
our effective model and the ab initio density functional 
calculation s "'^^i^^'^^i^^'^" can be seen from Fig. fUb) where 



we display the electronic dispersion from Ref. [15|. Note, 
the hole Fermi surfaces shifted by (tt, tt) is fully nested 
with that of the electron pockets which is also in full 
agreement with LDA results. Here, the position of the 
chemical potential fi has been deduced from the equa- 
tion n = A + x. We note that although the Fermi surface 
obtained previously in the effective two-band model^i re- 
produces correctly the one obtained within LDA calcu- 
lations, the actual evolution of the dispersion deviates 
significantly. 

Next we consider the one-loop contribution to the spin 
susceptibility that includes the intraband and the inter- 
band contributions: 




FIG. 1; (Color online) Calculated two-dimensional effective 
band structure along the main symmetry directions of the 
folded BZ (a) for the LaFeAsO system, and the reproduction 
(b) of the corresponding LDA band structurei^iiiJ^i^. The 
arrows in (a) and (b) indicate the points where bands cross 
the Fermi level. Note, the shade (color) of the curves in (b) is 
used just as a guide to the eye and does not reflect the actual 
hybridization of the bands. 




FIG. 2: (Color online) (a) Calculated Fermi surface topol- 
ogy for the LaFeAsO system. The arrows indicate the main 
scattering wave vectors. The filled dots refer to the states con- 
nected by the interband scattering at the AFM wave vector 
Qafm, while the open dots denote the interband and intra- 
band scattering at the incommensurate wave vector Qsdw- 
The dashed (cyan) lines and the +,— signs depict the po- 
sition of the nodes and the corresponding phase of the s± 
superconducting order parameter, (b) Calculated real part 
of the one-loop spin susceptibility along the main symmetry 
directions of the first folded BZ. The thick solid (black) curve 
refer to the total susceptibility while the other (red and blue) 
curves refer to the partial contributions which are multiplied 
by a factor of 5 for the sake of the presentation. The main 
scattering wave vectors shown in (a) are also indicated in (b) . 



and F'' are the normal and anomalous (superconducting) 
Green functions, respectively. 

In Fig. [2] we present the results for the real part 
of the total (physical) spin susceptibility xo(q;i'^m) = 
Si j Xo"* (li ^'^m)' ^c\\ as the partial contributions. 
The total susceptibility is dominated by the scattering 
at the AFM wave vector Qafm which is originated due 



Xo (q,i^m) 



2A^ 



to the interband {a — > /3) scattering. It is interesting to 
, i^ri)note that the intraband and interband scattering within 
a and /3 bands are very similar and are responsible for 
^2^he broad hump around the Qsdw wave vector. 

In the following we shall discuss the possible influence 
of the superconductivity driven by the short-range mag- 



+ F'{k + q, iw„ + iuJm)F^{k, iuJn)] 

where i, j again refer to the different band indices 
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nctic or charge fluctuations on the magnetic suscepti- 
bihty. It has been aheady argued that most Hkely the 
superconductivity in these family of compounds is of un- 
conventional origin and is driven either by the interband 
AFM fluctuations or by the intraband SDW fluctuations. 
However, one has to stress that even if the Cooper-pairing 
is driven by the interband fluctuations it still refers to 
the two fcrmionic states on the very same a or /3 bands. 
The standard Cooper-pairing for the two fermions from 
the different bands will be suppressed, since there are no 
states with k and — k that can be connected at the dif- 
ferent Fermi surfaces by the AFM momentum, as could 
easily be seen in Fig. [Ha). Therefore, we expect that 
inter-orbital AFM fluctuations will drive superconduc- 
tivity in the a and /3 bands. The latter should also re- 
sult in the very same value of the superconducting gap 
in both bands. The repulsive nature of the interaction 
would then rcquiro^^ the superconducting gap that sat- 
isfies At = — A{. , f-, . Thus, we consider the mag- 
netic susceptibility in the superconducting state assum- 
ing c?2.2_j^2-wave [Ak = ^ {coskx — cos ky)] and s±-wave 
[Ak = (coskx + cosfcj,)] symmetries of the order pa- 
rameter which both satisfy the condition given above. 

For the four-band model considered here the effective 
interaction will consist of the on-site Hubbard intraband 
repulsion U and the Hund's coupling J. There is also 
an interband Hubbard repulsion U' , which however does 
not contribute to the RPA susceptibility. Within RPA 
the spin response has a matrix form: 

Xi?,PA(q, iw,„) = [I - rxo(q, it^™)]"^ xo(q, i^^m) (3) 

where I is a unit matrix and xo(q, i'^^m) is 4 x 4 matrix 
formed by the interband and intraband bare susceptibil- 
ities determined by Eq. ^ . The vertex is given by 



non-SC, total Im^RPA 

■ cl^2_^z, total Im^RPA 

■ s±, total Im^RPA 



U J/2 J/2 J/2 

J/2 U J/2 J/2 

J/2 J/2 U J/2 

J/2 J/2 J/2 U 



(4) 



and we assume here J = Q.2U and U ^ t^^ . Note that the 
value of U was chosen in order to stay in the paramag- 
netic phase. We have to note that our interaction param- 
eters are carrying the band indices. Therefore, we neglect 
the possible orbital correlations. Whether this may play 
an important role needs to be addressed carefully and is 
beyond the scope of the present study. Though the cur- 
rent experimental and theoretical belief is such that the 
orbital physics is not involved in the physics of ferropnic- 
tides due to strong hybridization of all d-orbitals. 

In Fig. [3l[a) we show the results for the total RPA 
susceptibility, XflPA(q,it^rn) = Li j X^paI^' i'^™); as a 
function of frequency at the AFM momentum Qa_fa/- 
One finds that in the normal state the spin response 
does not show a well-defined peak but rather a broad 
continuum of the spin fluctuations. The origin for this 
is that the RPA enhancement of the AFM spin fluctu- 
ations is determined by the det |I — rx(q, itJm)|. One 





FIG. 3: (Color online) (a) Calculated imaginary part of the 
RPA spin susceptibility at the AFM wave vector Qafj\/ as 
a function of frequency in the normal and superconducting 
states. The thick dashed (red), dotted (blue), and solid (blue) 
curves correspond to the total RPA susceptibility. The thin 
(black) curves refer to the partial RPA contributions for the 
interband and intraband transitions in the s± superconduct- 
ing state, (b) Calculated imaginary part of the total RPA spin 
susceptibility in the s± state as a function of frequency and 
momentum along (1,1) direction. For the numerical purposes 
we set the damping constant 5^ — 0.8 meV. 



has to remember that the intraband on-site Coulomb re- 
pulsion U will strengthen the corresponding intraband 
fluctuations and the Hund's exchange will only increase 
directly the instability towards interorbital AFM fluc- 
tuations. Given the fact that each of the bare suscep- 
tibilities slightly differ from band to band as shown in 
Fig. [2tb), the RPA does not yield a well-dcflned pole. 
Thus one obtains simply a continuum of the fluctua- 
tions. The situation changes in the superconducting 
state. The quasiparticles at the Fermi surface of the a 
and /? bands connected by the AFM wave vector pos- 
sess the condition Ak = — Ak+q^^^jj- for the s± or- 
der parameter. The imaginary part of the interband 
magnetic susceptibility is zero for small frequencies due 
to the opening of the gap, and then it experiences a 
discontinuous jump at Q.c = min (|Ak| + jAk+q^^j^J) . 
Correspondingly, the real part of the interband [a — > 
(3) susceptibility will show the logarithmic singularity. 
This fulfils the resonance condition for the interband 
susceptibility: 1 — {J /2)Kcxq^ (SIafm ,^res) = and 
Itoxq^ {Q,AFM I'^res) = 0. Morcover, the intraband bare 
susceptibilities are small at this wave vector due to the 
direct gap, i.e. no states at the Fermi level can be 
connected by the Qafa/ for the intraband transitions. 
Therefore, a single resonant pole will occur for all com- 
ponents of the RPA spin susceptibility at Ures and 
the spin excition will form. This is evidently seen from 
Fig. [3l^a). Due to the single pole in the denominator 
all components of the RPA susceptibilities behave very 
similarly and the total susceptibility shows a well-defined 
resonance peak. 

In the case of dx2-y2 superconducting gap the situa- 
tion is more complicated. As clearly seen from Fig.l^l^a), 
the AFM wave vector connects states rather close to 
the node of the d^^^y^ superconducting order parame- 
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tcr and the overall gap in ImxQ determined by fic is 
small. At the same time even for this symmetry the 
resonance condition can be fulfilled due to the fact that 
Ak = — Ak+Q^j^j^. However, because of the smallness of 
^ Aq the total RPA susceptibility shows a moder- 
ate enhancement with respect to the normal state value, 
as seen in Fig. [3]Ja). Therefore, the resonance peak is 
pronounced only for the s± order parameter. Such a 
distinct behavior for the two various order parameters 
can be clearly resolved by the inelastic neutron scatter- 
ing experiments and therefore can be a direct tool to 
clarify the symmetry of the superconducting order pa- 
rameter in these systems. Like for d^2_y2 case, we have 
also found that there is no spin resonance for dxy- and 
dx2-y2 -l-idj^y-wave symmetries (due to their similarity we 
do not present these results). 

Finally we address the evolution of the resonance peak 
away from the AFM wave vector. In Fig. [3Jb) we show 
the total RPA susceptibility as a function of the momen- 
tum and frequency. Note that the s± superconducting 
gap changes only slightly at the a and /3 Fermi surfaces 
and can be considered nearly as a constant. Therefore, 
one always finds Ak = — Ak+q as long as the wave vec- 
tor q„ < Qafa/ connects the states at the Fermi surface 
of one of the a and one of the /3 bands. However, as 
it is also clearly seen from Fig. Eljb) the nesting condi- 
tion is very sensitive to the variation of q„ away from 
Qafm- Therefore, already at q„ « 0.995Qafm the 



Rcxq (q„, Wres) is much smaller than its value at Qafj\/- 
As a result the resonance peak is confined to the AFM 
momentum and does not disperse as it occurs for example 
in high-Tc cuprates. 

In conclusion, we have analyzed the behavior of the 
magnetic spin susceptibility in Fe-pnictide superconduc- 
tors. We show that the magnetic fluctuation spectrum 
calculated within RPA consists of (i) the continuum of 
the AFM spin fluctuations peaked at Qafm = (ti'iTt) 
that arise due to the interband scattering, and (ii) a low- 
q fluctuations around the C^sdw due to the intraband 
scattering. We show that the unconventional supercon- 
ductivity yields different feedback on the magnetic ex- 
citation spectrum. The resonance peak in form of the 
spin exciton occurs only for the interband scattering at 
the AFM momentum for the s± superconducting order 
parameter. We also find that the resonance peak is con- 
fined to the AFM wave vector and disappears rapidly for 
q < Qafj\/- 

Note added. After submission of our paper, we became 
aware of the experimental observation of the predicted 
resonance^^ and of the study by Meier and Scalapino^ 
who reached some similar conclusions as ours. 

We would like to thank A. Donkov, D. Parker, and 
P. Thalmeicr for useful discussions. I.E. acknowledges 
support from Volkswagen Foundation. 
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